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ABSTRACT 

The continuing development of the free electron laser (FEL) as a powerful and 
versatile source of coherent radiation steadily drives toward the goal of h.gh efficiency 
and broad tunability at shorter wavelengths. New experiments provide significant data 
and insight for analysis by theoreticians and experimentalists. Two important areas of 

study are short electron pulse effects, and the dynamics of optical mode distortion by 
intense electron beam currents. 

The initial part of this thesis examines one aspect of the projected task of FEL 
application as a military weapon. The advantages of the FEL over other directed 
energy sources are detailed, as well as the challenge presented by the effects of the 
marine atmosphere to high energy laser propagation. 

The remainder of this thesis examines several effects of long wavelength FELs. 
Chapter IV examines the proposed parameters of the CEBAF IR FEL, and the analysis 
leads to predictions describing system performance. Chapter V examines the effects 
of single pass optical mode distortion for FELs with narrow electron beams. Single- 
mode theory states that gain is proportional to the product of electron beam current 
and filling factor, but three dimensional simulations show that gain is a function of 
electron beam filling factor alone. Also examined is a phenomenon of destructive 
interference of light in the FEL undulator. 

Chapter VI extends the analysis of the relationship of gain and beam size to 
include multiple passes of laser light through the laser resonator. This affirms the 
general gain relationship, where gain is a function of electron beam filling factor, and 

also further explores the phenomenon of destructive interference within the optical 
mode. 
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I. INTRODUCTION 



A free electron laser (FEL) produces light through the stimulated emission 
from a high energy electron beam in a periodic magnetic field. Proposed by J. 
M. J. Madey in 1970 [1], the FEL has great potential for development as a 
highly efficient, powerful coherent light source with medical, industrial, and 
military applications. 

Theoretical models and technological demonstrations are extending the 
limits of FEL performance, with high power and efficiency possible within a few 
years. The FEL is a significant tool because of its unique operating 
characteristics. Among these are the use of electricity as the power source, 
order of magnitude tunability of output wavelength, and high "wallplug" 
efficiency. With the likely development of FEL power into the tens of megawatt 
range, the potential exists for military applications. Ground based and space 
based options have been studied by the Strategic Defense Initiative Office 
(SDIO). The option of arming a naval vessel with a high energy laser as an air 
defense weapon deserves study, and is discussed here in Chapter II. 

Chapter III provides a foundation of FEL theory to enable detailed 
descriptions of FEL performance in subsequent chapters. Chapter IV examines 
the infrared (IR) FEL proposed for the Continuous Electron Beam Accelerator 
Facility (CEBAF), of Newport News, Va. The planned accelerator and FEL 
parameters are analyzed to assist design performance optimization. The tuning 
range of the CEBAF IR FEL is analyzed to determine the influence of short 
electron pulse effects on FEL gain and power. 
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Chapter V treats optical diffraction, and the phenomenon of optical mode 
distortion caused by intense electron beams. Single-mode gain theory, which 
neglects the effects of diffraction and mode distortion, is shown to be 
inadequate. Using three-dimensional simulations of the optical mode wave 
front, a relationship between single-pass gain and electron beam size, 
accounting for diffraction and distortion, is found. A curve is generated which is 
generally applicable to all FELs with small electron beams. Chapter VI explores 
the new gain relationship and predicted destructive interference phenomena in 
a multiple-pass simulation. A similar gain curve is generated. 



2 



II. SHIPBOARD FREE ELECTRON LASER 



1. INTRODUCTION 

The current state-of-the-art anti-ship missile (ASM) threat posed to the U.S. 
Navy consists of high speed, maneuverable, robust and electronically 
sophisticated weapons. Current Fleet defensive weapons and tactics center on 
the highly complex Aegis Combat Systems Suite, which integrates sensors, fire 
control, and decision aids with surface-to-air missiles (SAMs). The Achilles’ 
Heel of this system is the limited number of SAMs and their performance 
limitations. High Energy Laser (HEL) design concepts have been studied for 
incorporation into combat ships for air defense. The Free Electron Laser (FEL), 
due to its design and performance characteristics, would be the optimal choice 
for a future ship based HEL weapon. 

2. THREAT ENVIRONMENT 

Combating the high speed and possible high density of anti-ship guided 
missiles is among the highest priority missions of tactical naval planning for the 
next generation of shipboard weapon systems. ASMs, capable of launch from 
submerged submarines, tactical aircraft, and small patrol boats are cheap, 
reliable and difficult to defeat with all but the most complex integrated systems. 
In worst-case wargaming scenarios, which include numerous ASMs 
simultaneously converging from many directions, even the most advanced air 
defense system, such as an Aegis-equipped guided missile cruiser, can be 
overwhelmed through exhaustion of the SAMs carried. This is known as the 
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"magazine depth.” 

Time is also considered a highly limiting factor for anti-ship missile 
defense. With present systems, once a target is identified, evaluated as a 
threat and the order is given to engage, there is a time lag due to SAM flight to 
intercept the target lasting from ten seconds for a close range engagement to 
two minutes for a longer range target. If the target is not destroyed, this 
process must be repeated, further wasting time and valuable defensive missiles. 
The Aegis system can successfully track and engage several incoming missiles 
simultaneously, but is ultimately limited by its magazine depth. 

3. HEL OPTIONS 

HELs, having been the subject of extensive research since the initiation of 
the Strategic Defense Initiative (SDI) program in 1983, are usually categorized 
by the fuel sources that drive them. These include X-Ray Lasers, Excimer or 
Chemical Lasers, and Free Electron Lasers. 

X-Ray Lasers, despite characteristics desirable for exoatmospheric SDI 
missions, are not a feasible option for a shipboard HEL. The wavelength of 
radiation produced by these lasers, in the X-Ray region, does not propagate 
effectively in the atmosphere, rendering it somewhat useless as a shipboard 
weapon. More importantly, the laser is powered by the detonation of a nuclear 
weapon, which understandably obviates further discussion of its adoption by the 
fleet. 

Chemical lasers, having been tested successfully as HELs, are a much 
more attractive option for a shipboard weapon. The MIRACL laser system, 
employing a deuterium fluoride (DF) laser, successfully destroyed a target 
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drone aircraft at White Sands Missile Range. This laser emits light in a 
favorable wavelength for atmospheric use ( 3.7-4.0n ) [2], and has been studied 
for possible shipboard use. The main obstacle for use of chemical lasers is the 
storage and combustive reaction of the highly volatile agents needed to drive 
the laser. The exhaust gases of chemical lasers are highly toxic and would 
require a complex air filtration and cleanup system to enable habitability for the 
crew of a vessel equipped with a chemical laser. 

Chemical lasers are driven through the exothermic reactions of highly 
volatile oxidizer and fuel streams which are mixed inside the optical cavity. The 
reactive nature which enables these agents to be useful is itself a strong reason 
for the rejection of chemical lasers as shipboard weapons. The transport and 
storage of highly pressurized, toxic, and reactive chemicals to ships at sea 
would require designs of replenishment ships (tankers and supply ships) that 
today do not exist. The potential for disaster once the chemicals were aboard 
the laser ship would necessitate detailed and expensive safety programs, 
further complicating the process. Ultimately, the chemically driven laser system 
would fall victim to the same limitations of the missile systems it would seek to 
replace. The chemical laser is driven by an irreversible chemical reaction that 
would cease once the fuel and oxidizer were exhausted. This is a solid 
magazine depth that would severely limit the duration of engagement and 
require complex reloading procedures. This is a major flaw for a plan to field a 
chemical HEL aboard a ship. A principle advantage of a ship as a weapons 
platform includes long endurance on station, and it would be wasteful of this 
endurance potential to place a short-legged weapon aboard as a principle 
means of air defense for a battle group or geographic region. 
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A simpler and more logical system for deployment of HELs for near and 
long term use would involve development of an aircraft-based chemical HEL, 
followed later by the development of a shipboard FEL (SFEL) when the 
necessary technology matures. Aircraft lack the endurance of ships, needing to 
land frequently for fuel, maintenance and crew rest, so that the engage-and- 
reload requirements for a chemical HEL would not represent as severe a 
mission limiting factor as would be the case for a shipboard chemical HEL The 
lessons learned through development of the aircraft based HEL would later 
benefit the deployment of the SFEL 

The Free Electron Laser has several advantages over other potential HEL 
designs. The principle advantage of this system is the wide range of 
wavelengths over which the FEL can be tuned. This allows the FEL to be 
tuned to exploit highly variable atmospheric conditions and to be more highly 
resistant to countermeasures. An additional advantage, crucial when 
considering distribution of large quantities of electricity aboard ships, is that 
scaling up to allow the addition of an FEL would be simply a matter of adding 
more generators and power conditioning equipment. Recent technological 
advances in capacitor storage banks would allow large quantities of electrical 
energy to be held in reserve for immediate use, providing a quick reaction 
capability for the SFEL. This is an easily manageable task, and would fall in 
well with the navy’s future plans for an integrated electric drive combatant 
vessel. An especially attractive factor of FEL design is the generally accepted 
"wallplug" efficiency goal of 40%. This is would be made possible through use 
of an energy recovery system, which would also enable reduction of the 
necessary mass of radiation shielding. The limiting factor for the operation of 
the FEL would then be the common fuel supply used for propulsion and 
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electricity of the ship, which is already supported by a well established system. 



4. ATMOSPHERIC PROPAGATION EFFECTS 

The marine atmospheric environment within which an SFEL would operate 
is a highly variable, dynamic system with numerous linear and nonlinear effects 
that have profound influence on the propagation of high-power laser beams. In 
addition to the effect of optical diffraction within the atmosphere, there exist 
atmospheric attenuation effects. The mechanisms of laser beam photon loss 
through atmospheric attenuation by the constituents suspended in the 
atmosphere can be broken down into several distinct effects: the linear 
absorption and scattering of the molecular and aerosol elements in the 
atmosphere, the random wander, spreading, and distortion of the laser beam 
due to the effects of turbulence within the atmosphere, and the non-linear 
defocusing effects of thermal blooming caused by the absorption of small 
amounts of beam power by the atmosphere. These effects are dependent on 
atmospheric temperature, clarity, humidity, wind velocity, and beam 
characteristics such as wavelength, temporal operating mode and beam 
diameter. 

Accounting solely for linear atmospheric effects, the mean peak irradiance 
of a propagating laser beam is l P = [P/{nw 2 )] exp (-a r z), where P is the 
transmitted laser power, z is the range, and a r is the atmospheric extinction or 
attenuation coefficient. The term w is the root mean square (rms) beam radius 
and is a function of many factors, w = yfwj + w 2 + w f 2 , where w d is a factor of 
diffraction and beam quality, w y accounts for beam jitter, and w, is a correction 
for atmospheric turbulence [3]. 
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The diffraction and beam quality factor correction is 
w d = Vp 2 z 2 //c 2 w 2 + w 2 ( 1 - z/R) 2 . Here p is a ratio of actual focal spot radius to 
the diffraction limited spot radius, or 



o _ 1/e Radius of actual t££§i spot _ w d _ . « 

P 1/e Radius of Diff. Limited spot w /dea/ ‘ ,d9a/ 



2-1 



Range R can be varied from R = z for a beam focussed to the point at 
distance z, to R = ~ for a collimated beam. Beam jitter effects are corrected for 
by the addition of the mean square radial displacement of the focal spot, 
Wj = V2<9 J>P, at range z for a varying single axis jitter angle, <e£>. This angle 
is isotropic for x and y axes, so that <9|> = <9 2 >. 

The attenuation and extinction term, a T = a abs + a scat , represents the 
effects of absorption and scattering of molecular and aerosol components of the 
atmosphere. Scattering is a directional redistribution of the beam incident on 
the various molecules composing the atmosphere. This does not reduce the 
energy of the beam but does reduce the intensity. Absorption is a very 
complicated and irregular function of beam wavelength, and is a critical factor in 
the effects of thermal blooming. Scattering and absorption by aerosol 
components is generally less dependent on wavelength than molecular 
absorption and scattering. The slight wavelength dependence of the absorption 
and scattering coefficients is a result of the combined effects of the aerosol 
complex refractive index and the size distribution of the aerosol particles. 
Aerosol effects become somewhat significant at shorter wavelengths 
(X<3|4. m)[3j. Fig 2-1 shows the highly variable nature of the aerosol and 
molecular effects. 
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Figure 2-1: Calculated aerosol absorption and extinction by the 
atmosphere in 23km visibility vs. wavelength [3]. 



Figure 2-1 is generated through a model of the continental atmosphere. 
The marine environment would include higher levels of humidity and its effects. 
Included on the graph are range bars indicating the wavelengths of three 
common lasers and the variation of the attenuation coefficients for each. Of 
note is the dip in aerosol absorption effects for A. = 3.8nm, which closely 
corresponds to the wavelength of the DF laser. This wavelength is near an 
optimum wavelength selected for the study of FEL placement on a naval vessel 
[4]. 

Atmospheric transmittance, accounting for molecular and aerosol 
absorption and scattering effects, is a ratio of initial irradiance, / 0 , to that at a 
range z : 
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Atmospheric turbulence arises from convective air currents and eddies 
caused by air-surface temperature gradients, which result in slight fluctuations 
in the index of refraction. For beam propagation, the eddies form pockets of 
disturbed air, known as "turbules", which vary in size relative to the beam 
diameter. This non-uniformity results in the various effects of beam wandering 
off the target, intensity fluctuations across the beam profile, and spreading of 
the beam spot size. Turbulent beam spread and wander effects are accounted 
for by the term w, = 2zikL 0 where z is distance, k is the wave number, and 
L 0 = (0.545k 2 C^z) _3/5 is the lateral coherence length of a spherical wave. Cft is 
the refractive index structure constant. Figure 2-2 shows wavelength 
dependence of the effects of beam spreading and wander on average 
irradiance, l(k), normalized by the irradiance at X = lOja.m [3]. At long ranges, 
turbulent beam spreading becomes wavelength independent, but this is an 
exceptional regime. In the usual case, the turbulence coefficient 
w, = 2.0l(r 1/5 c£' 5 z 8/5 ). Though the dependence of turbulent beam spread on 
X -1 ' 5 is relatively weak, this relationship implies that there is an optimum 
wavelength for propagation for given values of turbulence level, range, and 
aperture diameter (diffraction being proportional to X). Neglecting jitter effects, 
the mean peak intensity is given by 



1 1 
w d 1 + W t 2 lw$ 



2-3 



Assuming a reference wavelength Aq, the ratio of the irradiances 

/(X) = f^o! 2 1 + A 

l( Ao) LM 1 + A(AcA) 12/5 ' 



where A = 16.l2jr 2 w 2 (C^zA < j 2 )'* /5 /p 2 . For wavelengths shorter than the critical 
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wavelength, defined by X c = (/A/S) 5 ' 12 ^, the effects of turbulent spreading 
dominate the diffraction limited focusing, and beam intensity decreases as X 2 * 5 
[3]. 




Figure 2-2: Dependence of turbulence effects on wavelength [3]. 



Thermal blooming or defocusing is a time dependent effect of absorption of 
a small amount of beam energy by the atmosphere along the beam path. The 
resultant rapid heating lowers the air density and also the refractive index, 
causing a nonlinear negative lens effect of divergence on the beam. Varying 
with atmospheric conditions, defocusing can occur in as little as 0.2s. This 
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length of time is crucial because it is close to a typical value of the "dwell time" 
necessary for a laser beam to remain focussed in order to damage or destroy a 
target. The spreading effect can dominate diffraction limited focusing. This can 
limit the irradiance on a target independent of the available power of the beam. 
Accounting for variable atmospheric velocity due to wind and the slew rate of a 
beam tracking a target, there will be a temperature gradient across the face of 
the beam, with the "cooler" air on the "upwind" side of the beam face, in the 
direction of the cumulative wind/slew velocity. The temperature gradient will 
result in a varying refractive index across the beam face, leading to an 
asymmetric distortion and growth of the beam cross section. With the "upwind” 
side of the beam face remaining at a lower average temperature, this side will 
tend to have a higher air density and thus a higher index of refraction. This 
causes the path of the defocused beam to bend "upwind” asymmetrically. 

Peak irradiance on a target accounting for thermal bloom effects must be 
scaled by a factor of thermal distortion, N, which is a function of focusing of 
beam radii, the effects of a non-uniform velocity profile, and atmospheric 
attenuation, a T . The specific heat, c p , of the air, the range, z, and the power, 
P, of the laser are also accounted for. Essentially, N « a T PoZ/vD 3 , where v is 
a wind velocity term, a 0 is the collimated beam 1/e radius, and D = 2 V 2 a 0 is the 
aperture diameter. An empirical relationship that works well in modeling 
systems evolves as l rei = Ibioom/lpeak = 1/(1 + 0.06/V 2 ). Since peak irradiance is 
proportional to the product of power and the empirical correction, l peak « Pxl rgh 
the peak irradiance on target will increase until l rei decreases at a rate 
proportional to 1/P or greater. This indicates that for a given beam size, focal 
range and atmospheric conditions, there exists a peak power, P c , at which 
irradiance on target is maximum. 
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Using this empirical relationship, Fig 2-3 shows the effects of thermal 
blooming in limiting peak irradiance and maximum power for a propagating 
beam. 




P/Pc 



Figure 2-3: The normalized peak irradiance vs. the normalized power PiP c , 
based on the empirical relationship of Ib/oomUpeak [3]. 

Unlike other atmospheric effects, increased power in the laser beam will 

not result in higher fluence on the target, but will only serve to further aggravate 

the blooming effect. This is a nonlinearity unlike any other atmospheric effect. 

All other phenomena result in increased target irradiation with higher power. 

Peak irradiance accounting for linear and nonlinear effects at range R = z is 

P expHtrz) 2-5 

na 2 1 + 0.0625 N 2 
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The linear effects of diffraction, beam quality, jitter, and turbulence are 
accounted for by the rms beam radius 

w 2 = 2 P 2x2 f 2 + 2<e 2 >z 2 + 4(c^z) 6/5 r 2/5 z 2 2-6 

n 2 D 2 

where D = 2V2a 0 is the aperture diameter, and the linear attenuation effects are 
covered by a r = a abs + a scat . 

The effects of thermal blooming are time dependent, and so there is a 
difference between the propagation effects for continuous wave (cw) and pulsed 
beams. For a pulsed beam of diameter D, if the pulse repetition frequency 
(PRF) is less than the time needed for the cross wind, v 0 , to clear the tube of 
atmosphere heated by the previous pulse, or PRF < v 0 id, there will be no 
thermal bloom effects from overlapping pulse beam effects. For time-scales 
outside of this regime, however, both pulsed and cw beams will experience 
bloom effects due to continuous heating. Through complex studies of various 
PRFs and beam face irradiance levels for pulsed beams, it has been found that 
pulsed beams are less affected by atmospherics than cw beams at long ranges 
[3]. 

Marine atmospheric propagation effects are difficult to predict on any useful 
time scale. The few controllable parameters include light wavelength and 
power level of the laser. Chemical HEL designs studied for naval application 
have all had basic flaws of limited fuel capacity and fixed wavelength. Free 
electron lasers, though less mature in design, have the valuable characteristics 
of tunable wavelength and manageable fuel and power requirements. The 
minimizing of atmospheric effects on laser beam propagation will need to be 
carried out not only by the selection of laser parameters, but also in creative 
beam optics. A system capable of atmospheric propagation diagnostics and 
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adaptation of deformable optics would greatly enhance the performance of any 
high-energy-laser weapon. Atmospheric diagnostics, hitherto the realm of 
astronomical studies for adaptive telescope optics, should be developed for 
inclusion into the control system of any SHEL design. In astronomy, optical 
system diagnostics are performed by observation of a ’’guide star" of known 
characteristics. For tactical military purposes, a low power laser carried by 
drone aircraft or other suitable platform could be aimed to illuminate diagnostic 
optics to allow analysis of the prevailing atmospheric conditions and subsequent 
optimization of the beam parameters and adaptive optics’ adjustment [2,5]. 

A more exotic solution may be found by stationing a relay mirror system 
high above the FEL vessel aboard a long-endurance fixed-wing remotely-piloted 
vehicle (RPV) or balloon aircraft. A system centered on this concept, named 
Thunderball, was presented to the Physics department of the Naval 
Postgraduate School by Lieutenant Colonel Ed Pogue, USA, former Deputy 
Director of Directed Energy for the Strategic Defense Initiative Office. This 
concept had been studied by the Advanced Technology Group of W. J. Schafer 
Associates [6]. The chief advantages of this system, whose aircraft would be 
stationed at altitudes approaching eighty thousand feet, would be the 
minimization of distance traveled through the dense lower atmosphere. This 
would help minimize the most harmful effects of the atmosphere. Additionally a 
beam of greater diameter, and thus lower power density across the beam face, 
could be used. This would lessen the damaging effects of atmospheric 
influences. Once the beam reached the optical system aboard the aircraft, a 
beam cleanup and refocusing system could narrow and redirect the beam 
toward a target far beyond the horizon of the FEL vessel below. This system 
would be well suited for a defensive mission against incoming ballistic missiles. 
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III. FREE ELECTRON LASER THEORY 



A. BASIC FREE ELECTRON LASER PHYSICS 

Free electron laser systems consist of four basic components, the electron 
beam accelerator, the optics which include the beam transport system and the 
FEL undulator and resonator. The accelerator produces a stream of relativistic 
electrons that is injected into the periodic magnetic field of the FEL undulator, 
which "wiggles" the electrons along sinusoidal trajectories. The periodic 
transverse acceleration of the highly relativistic electrons results in classical 
synchrotron radiation into a forward cone along the beam path. The resonator 
cavity consists of two curved mirrors, one is fully reflecting and one is partially 
reflecting, placed beyond the ends of the axis of the undulator. The mirrors are 
designed to sustain a transverse Gaussian mode shape. The total resonator 
length is typically S = 5 - 20m. The FEL oscillator starts from spontaneous 
emission, some of which is saved in the resonating cavity resulting in an 
increase in optical power. Repeated coupling of the electrons with the optical 
and undulator fields results in the bunching of the electrons. This amplifies the 
amplitude of the wavefront, resulting in coherent radiation and lasing. Undulator 
gain is the fractional change in optical power per pass, G = A P/P. Typical 
values for G = 0.10 for a low gain FEL and G = 10 4 for a high gain FEL. 

The bunching effect is crucial to FEL operation, and is highly dependent on 
the quality of the electron beam produced by the electron accelerator. 
Specifically, the energy spread of the electrons as well as the angular spread of 
the beam must be reasonably low for FEL success, the energy spread among 
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the electrons in the beam results in differences of velocity within the undulator, 
leading to an axial smearing of the electron bunch entering the undulator. If the 
velocity differences spreads the electrons beyond the optical wavelength, X, the 
FEL gain would be seriously degraded. Variations of the electron injection 
angle, or angular spread, results in non-axial components of bunch velocity. 
Reducing the axial velocity will also serve to degrade gain in much the same 
way as an energy spread. 

Emittance and beam spreading, collectively known as "beam quality," 
influence bunching effects. Emittance, e = re, is the product of the root-mean- 
square (rms) values of radial and angular spreads of the electron beam. In the 
tuning of a given electron transport system, both of these parameters may be 
changed, but their product, e = r0, remains constant. The normalized 
emittance, e„ = ye, is the product of the electron beam Lorentz factor, y, and the 
emittance e. The normalized emittance tends to remain constant even when 
the electron beam is accelerated to higher energies. Normalized emittances 
typically vary from a few tens to a hundred mm-mrad. 

The accelerator produces relativistic bunches of electrons with peak 
currents from 1 A to IkA, and bunch lengths of l e = I0 _1 cm. The relativistic 
electron beam will typically have a Lorentz factor, y, from ten to the low 
hundreds. Typical electron beam radii are r g = 0.10cm. These result in electron 
densities of p = lO^cm -3 to 10 13 cm -3 . With Lorentz factors y varying from ten a 
few hundred, the beam energies, ( y - 1 )mc z , range from several tens of to 
hundreds of MeV, where m is the electron mass, c is the speed of light, and 
me 2 = 0.51 1 MeV is the electron rest mass. 
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The undulator consists of a series of alternating magnetic pole pieces 
forming an intense magnetic field that accelerates the relativistic electron beam 
into periodic transverse motion. Within the undulator, electron bunches undergo 
one cycle along an axial length known as an undulator wavelength, Xo, which is 
typically five to ten cm. The number of undulator wavelengths, N, varies from 
twenty five to 250, which results in undulator length L = NX 0 from one to ten 
meters. The undulator magnetic field strength, S, ranges from 3kG to 5kG. 
Mean undulator magnetic field strength, B, is rms field averaged over the 
undulator period X^. For helically (circularly) polarized undulators, B = B, due to 
the constant acceleration of the electrons in a helical path along the axis of the 
undulator. Helical undulators produce circularly-polarized radiation. Magnetic 
fields for linearly-polarized undulators result in non-constant acceleration due to 
the sinusoidal electron path, so that B = B/V2. Linearly-polarized undulators 
produce linearly-polarized radiation. 

These parameters, with the electron charge magnitude, e, combine as the 
"undulator parameter" [7], 



The undulator parameter, K varies from 0.1 to 3. K may be varied by altering 
the separation distance between pole faces along the undulator in order to 
change B. The optical radiation wavelength, X, produced by the FEL is given 
by [7] 



The optical wavelength is a function of the undulator magnetic field and the 
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electron energy from the accelerator, as represented by the Lorentz factor y. 
This relation reflects the tunability of the FEL over a very broad range of 
wavelengths. Experiments have a xIO range of tunability, and the FEL can 
operate from the mm to the X-Ray wavelength range. 

Short optical wavelengths are a result of the Doppler shift effect on the 
undulator and optical radiation fields as "seen" by the highly relativistic 
electrons. The electrons "see" a Lorentz contracted undulator wavelength, 
Xo ~ Xo /y. The electrons emit light at a longer Doppler shifted wavelength 
X' = 2yX. With the similarity of the Doppler shifted light and the Lorentz 
contracted undulator wavelength, X' = V. the relationship of the radiation 
wavelength to the undulator wavelength becomes X = Xo/2-y 2 . This is the 
mechanism which enables short wavelengths to be produced. 

The resonance condition is met when one wavelength of light passes over 
an electron as the electron traverses through the length of one undulator 
wavelength. As the electron bunch is subjected to the resonant forces of the 
combined undulator and optical fields, the electrons undergo either stimulated 
emission or absorption, depending on the phase relationship of the bunch 
relative to the optical field. This new light from the electrons will interfere with 
the light of the optical mode. If the net interaction of the radiation is 
constructive interference, there exists a regime of positive gain, G = a P/P > 0, 
for the undulator. 

Figure 3-1 illustrates the interaction of an electron with the optical mode 
and magnetic fields within the undulator. The top diagram shows an electron 
bunch and an optical pulse in the undulator. The middle diagram illustrates the 
details along one undulator period. 
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Figure 3-1: Fundamental design and electron-optical mode interactions 
within the undulator of an FEL. 
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The electrons follow a sinusoidal trajectory along the periodic magnetic 
undulator field while one wavelength of light passes over each electron. The 
bottom diagram details the electron-optical mode interaction. Depending on the 
relative phase between the optical mode and the electron, the electron will be 
subject to either a retarding or accelerating force that will induce it to either emit 
or absorb energy. Conservation of energy determines that the energy 
transferred from the electrons must be absorbed by the optical mode wave. 
This results in optical amplitude growth and FEL gain. Should the phase 
relationship cause the electron to be accelerated, the force to do so will rob the 
optical mode of energy and result in decreased amplitude, or negative FEL 
gain. In non-resonant systems, gain can be positive or negative in various 
positions along the undulator resulting in complex optical amplitude profiles. 



B. ELECTRON DYNAMICS 

The highly relativistic electron beam is subject to the combined forces of 
the radiation of the optical mode and the magnetic field of the undulator. The 
equations of the Lorentz force exerted on the electrons are 

? = = -e [£ + $ x £f)j, (y me 2 ) = -ec(j? • £) 3-1 

where $ = V/c is the ratio of electron velocity V to the speed of light c, £ is the 
optical electric field, and the net magnetic field, £ = £ r + £ 0 is the sum of the 
magnetic fields of the optical mode and the helical undulator [7]. Using the 
combined undulator and optical fields in 3-1, and integrating over time, the 
transverse electron motion within the undulator may be solved. There are a 
total of five equations in 3-1, with only four unknown terms p?(f). 7(0), so that 
one of the five equations may be disregarded. The relativistic electrons have 
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an initial position z 0 , an axial velocity V z =$ z c, where p z = 1 - (i + K 2 )i-f = 1, 

and the position along the undulator axis is given by z{t) ~ z 0 + $ 2 ct + The 

transverse velocity p a = Vi - Pz « 1 . so the magnitude of the transverse motion 
of the electrons is much less than the axial motion. The transverse undulator 
fields are much greater than the transverse optical fields, so that the transverse 
force equation of the electrons becomes 

?1 = “tCy Pi) = -=Jpz S[-sin(/c 0 z), cos{k 0 z), 0]. 3-2 

at me 

Where k 0 = 2tVXq is the undulator wave number. Integrating this equation over 
time, and assuming perfect beam injection, yields the transverse velocity: 

ft. = ^ [ cos {k 0 z), sin(/c 0 z), 0]. 3-3 

One more integration will result in the transverse position relationship for the 
electrons: 

= = 1 ^ [ sHk ° z) ’ ~ C0S{k ° z) ' 0] - 

It can easily be seen here that because K = 1 and the Lorentz factor y= 100 
that the magnitude of the transverse motion is relatively small. 

C. ENERGY EVOLUTION 

The second term in the Lorentz force equations 3-1 determines the 
evolution of the energy. Using the expression above for transverse velocity p 4 
3-3, the approximate electron position z = z 0 + po ct + • • • , for Lorentz factor 
1, the energy for a weak field regime may be approximately calculated. The 
optical electric field is 
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[z r = E (cos(kz - tof), -sin (kz - of), 0), 

where the wave number is given by k = 2nfk, and the principle carrier frequency 
is given by (o = kc = 2nd\\\. Using y = -efl? • lz)/mc, the initial position 
z 0 = z 0 + p 0 cf, ar| d substituting the expression for £ r , the energy term evolves 
as 

y = — cos(( k + k 0 ) z - (of) = — — cos(( k + k 0 ) z 0 + [( k + k 0 ) Pq - k] ct), 3-4 

ymc y me 

where p 0 is the ratio of initial electron velocity v 0 and the speed of light c [8]. 
Define the initial phase as Co = (^o + k ) z o = 2jt z o / ' k > which relates the initial 
electon z position relative to the optical wave. The time dependent phase is 
C(x) = {k 0 + k)z{t) - cof. Electrons separated by AC, = 2n are in identical positions 
in adjacent optical wavelengths. The phase velocity is v 0 = L[{ k 0 + /()p 0 - k], 
and the dimensionless time is t = cf/L. The energy evolution may now be 
written as y = {ekE/ymc) cos(£ 0 + v 0 t). This can also be restated in a 
dimensionless time regime where (.°.) represents a dimensionless time 

o 

derivative d(...)/dx, so that y = {eEKL /ymc 2 ) cos(Co + v 0 t). Integration of this term 
will result in 

Y = Yo + — [sin(£ 0 + v 0 x) - sin(^o)], 

Yo mcS'o 

where y 0 is the initial energy. Rearranging, the time varying energy in the weak 
optical field regime is 

Y(t) = Yo + — [ sin(Co + v 0 t) - sin(Co) ] , 
v 0 

where the term A = eKEL/y 0 me 2 determines the energy exchange. This shows 
that the energy exchange of an FEL is controlled by the initial energy Yo of the 
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electrons, the magnitude of the electric field E, the magnetic field strength 
K « b, ant the electron phase t^ 0 = 2nz 0 /\ within each optical wavelength, and 
the electron phase velocity v 0 . For large v 0 (vo>*7c) the change in electron energy 
oscillates rapidly along the undulator and is small. For v 0 near resonance 
(v 0 = 0), the electron energy changes are linear in x, y(t) = y 0 + Axcos(Co). and 
are periodic in the electron phase Co- For -7t/2<Co<rc/2, the electron energy 
increases, and for 7t/2<Co<3rc/2, the electron energy decreases. For the whole 
beam randomly spread in Co. the average energy change is zero. To this order 
in the field E, there is no gain. 



D. THE PENDULUM EQUATION 



The FEL phase and phase velocity relationship may be calculated by 
returning to the Lorentz energy expression, y=-{e/mc) • £ f . FEL gain is the 

net energy exchange to the optical field and must be calculated to higher order 
in the field E. Substitution of the optical mode electric field into 3-4 results in 
an expression of energy gain [8] 



1 = 



y 



eKEL 

■fmc 



cos(C + <t>). 



To track the electron phase changes as the energy changes, the Lorentz factor 
is rewritten y 2 = 1 - f 2 = 1 - (pf + p|) where Pj. = K/y. For highly relativistic 
electrons y » 1 , and p z = 1 - (1 + K 2 )y 2 . Phase velocity t, = v = (k 0 + k)z - co, 
and the second derivative of phase, the rate of phase velocity change, is 

C = v = (* 0 + k)z = kc = kc cos(C + 4>). 3-5 



Equation 3-5 has the form of the Pendulum equation. Taking advantage of the 
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fact that the initial carrier frequency too = k 0 c = co (1 + K 2 )/2-f, 3-5 can be 
restated as 



2eKEk 0 

■fm 



cos (£ + <{)). 



The dimensionless time x= ct/L, wherein 0< x < 1 as the electrons travel down 
the axis of the undulator, may be inserted into these expressions. This requires 
that d(..)/dt = L/c d(..)/di or (..) becomes (°.), and 

oo 

d 2 (..)/dt 2 = (L/c) 2 d 2 (..)/dt 2 = (..). This now allows the Pendulum equation, 3-5 
be rewritten: 



°o 0 I . 1 ^ 2.&EKk 

C(x) = v=^ — [• — 7 — - cos(C + <M = la I cos(C + <{>), 3-6 

l c J 

where the dimensionless term la I = AnNeKEL/y^mc 2 is the strength or 
magnitude of the optical field for a circularly-polarized undulator [9]. la I shows 
that higher velocity (larger y) electrons would require a stronger optical field E in 
order to achieve the same electron bunching. When lalcrc the field is 
considered to be in a "weak field" regime, and the electron bunching after a 
single pass through the undulator is negligible. When la I = 7x the bunching 
effect becomes noticeable. While la l=»>7x is the "strong field" case, where the 
electrons quickly bunch within the undulator and may cause the onset of 
"trapped particle instability." If the optical power is strong enough to cause one 
synchrotron oscillation for each pass, the electrons become trapped in deep 
potential wells. This results in beam current oscillation that destabilizes the 
carrier wave, causing sideband frequencies to grow. The coherent wave that 
had been established by the FEL interaction, may be amplified or degraded by 
the oscillating current of the trapped electrons. [9] 
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E. THE WAVE EQUATION 



The FEL oscillator starts with spontaneous emission and creates a 
classical light wave that is amplified by stimulated emission. Because the 
spontaneous emission bandwidth of the wave is approximately 1/A/, where N is 
the number of periods in the undulator, there is some measure of coherence in 
the initial light wave. Since A/»1 , the bandwidth is sufficiently narrow to form a 
classical wave with a complex field envelope (E(z,/),<j>(z,0) whose circularly- 
polarized electric and magnetic fields are expressed by 

= E(z,f)[cos(vy) ,-sin(\y),0], #,(*.0 = E(z, f)[sin(\y),cos(\|/),0] 3-7 

where E{z,t) is the wave amplitude, y = kz - cot + <t»(z ,f), and the phase is <t>(z,f). 
The amplitude and phase vary slowly in time over the optical period 
(E«ccoE,<t»-cax}>), and also in space over the length of an optical wavelength X, 
such that (EckE,( K/«t>) where k = 2kJ\ [8]. This approximation allows the fields 
to be derived from the vector potential [9]: 

/?(*./) = -^^[sin(v),cos(y),0]. 



The optical wave equation becomes 



= 2 



dE J dE 
dz c dt 



[cos(\j/),-sin(v|/),0] + 2 E 



3<j> 1 3<t> 

dz c dt 



[-sin(\|/),-cos(v|/),0] 



= ^7, (*,/), 3-8 

c 

where d 2 = (P 2 - (He 2 ) (d 2 /dt 2 )) is the D’Alembertian operator, and is the 
transverse beam current. If there is no (no transverse motion of the 
electrons), then the solution to 3-7 is E = E(z-ct) and <|> = <t>(z-cf). 
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Each electron contributes to the total transverse current so that 



7 X = -eclj?j8 (3) (x* - T,), 



3-9 



where 7- is the trajectory of the i* electron and 8 (3) is a 3 dimensional Dirac 
Delta Function [8]. The coherence volume, V c , is the product of the gain 
bandwidth length, NX, and the transverse mode area, XL=XNXq, so that 
V c = N 2 X 2 Xq. In order to correct for the fast variation of the trigonometric 
functions of y in 3-7, an average is taken over a small volume element dV c 
which is small relative to the coherence volume, but large relative to an optical 
wavelength X. [9] Using this averaging and inserting 3-7, 3-8 becomes 



JL + ±i_ 

dz C dt 



Ee '* = -4neKp(z - p z ct)<e '^> 



3-10 



where p(z - (? z cf) is the mean density of the electron bunch. If the longitudinal 
spread of the bunch is so great that the density is no longer a function of axial 
position, there will be uniform evolution of the electron pulse along its length, 
and no optical modes are followed, so that 3-10 is no longer dependent on z. If 
the number of undulator periods N is large, there is a small bandwidth and y 
changes only slightly, and the wave equation becomes 

a = -j<e-^> 3-1 1 



where 



4k NeKELe^ . _ 8/V(e;tKL) 2 )p 
■fmc 2 ’ 1 y*mc 2 

The field magnitude is a, and j is the dimensionless beam current [8,9]. The 
beam current, j, determines the reaction of the optical mode to the bunching of 
the electron beam, and provides coupling between the wave and the electron 
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beam. When j < n the gain is low, and when j the gain is high. 

For a linearly-polarized undulator, the non-uniformity of electron 
acceleration requires that these expressions be modified by K^>K(J 0 {ty - J^)). 
Here K is the undulator parameter, £ = K 2 / 2(1 + K 2 ), and J 0 and J ^ are Bessel 
functions of the first and second kind. [9] 

F. DIMENSIONLESS PARAMETERS 

Certain dimensionless parameters recur throughout the descriptions of the 
physics of the FEL. Their use allows quick summarization of performance 
without regression to involved calculations. Several important dimensionless 
parameters have already been defined and are used in this paper. They 
include the undulator parameter, K, the electron phase and phase velocity 
C and v, the beam current j, and the optical field magnitude la I. In the 
derivation of 3-10, it has been assumed that electron and optical mode sizes 
are the same. The filling factor F = n,r 2 /nW$ is the ratio of the area of the 
electron beam to that of the optical mode waist. When the electron beam is 
entirely within the optical mode, the dimensionless beam current j (3-10) should 
also include the filling factor F as a correction for the smaller volume of the 
electron beam. 

The transverse dimensions can be normalized by the optical mode radius 
v'rt/L X. Macroscopic longitudinal dimensions can be normalized by the undulator 
length L. The normalized electron pulse length is c z = l e iN\ where l e is the 
electron bunch length. The Rayleigh length Z 0 is a measure of optical 
diffraction properties and is dependent on optical wavelength, mirror design and 
configuration. The normalized Rayleigh length is z 0 = Zq/L. Z 0 is related to the 
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optical mode waist radius, W 0 , by Z 0 = kWq/X. The normalized mode waist 
radius, w 0 = W^kL/X, is related to the normalized Rayleigh length by 
w 0 = Wq'IkL/X = 'JXZq/k Vj iL/X = JT 0 . 

The slippage distance As is the difference in distances traveled by the 
optical wave and the electron beam as they move down the undulator. This is 
also the length that the light moving at c passes over the slower electrons 
moving at (3 z c during this period. Their relative speed, Av = (1 - p z )c, multiplied 
by the period of one undulator pass At = Lie, is the slippage distance 
as = Av At = (1 - p z )c(L/c) = (1 - p z )L. Inserting the definition of p z , undulator 
length L = NXq, and wavelength X = Xq(1 + K 2 )i2-f, the expression for slippage 
distance becomes as = NX. The normalized electron pulse length c z = l e /NX is 
scaled by the slippage distance [9], 
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IV. THE CEBAF INFRARED FEL DESIGN 



A. BACKGROUND 

The Continuous Electron Beam Accelerator Facility (CEBAF) in Newport 
News Va. is administered by the Southeastern Universities Research 
Association (SURA) for application to nuclear science, technology and industry. 
The CEBAF superconducting accelerator incorporates two 400 MeV linear 
accelerators (linacs) with interconnecting, recirculating electron beam lines. 
These will allow five simultaneous electron beams at energies between 800 
MeV and 4 GeV to be delivered to three end-stations. A program is underway 
to design and construct two high power FELs operating at the near infrared (IR) 
and the deep ultraviolate (UV) wavelengths. The design configuration of the 
FEL facility will allow laser operation without interference to the nuclear science 
research. The IR FEL will take advantage of the operational 45 MeV 
accelerator injector as its driver, while the UV FEL will use the 400 MeV North 
Linac. [10] 

The CEBAF facility has several features which make it interesting for 
FELs. The CEBAF linacs will supply the FELs with a high-quality electron 
beam with normalized emittance e„ = 15 mm-mrad, a high average power of 
800 kW, energy spread Ay/y = .001, and a duty cycle of 100 %. The fraction of 
energy extracted by the FEL is 1 1N, so that the average power of the FEL 
would be = 30 kW. 
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Figure 4-1: CEBAF IR and UV FEL configuration. The UV FEL will 
utilize the 400 MeV North Linac for its driver, while the IR FEL will use 
the 45 MeV Injector [10]. 



These are important parameters because if the FEL were ever to be developed 
for military use, high average power and performance goals would have to be 
achieved, and this experiment will be an important step. Additionally, the 
accelerator uses advanced superconducting technology that will be incorporated 
into the next generation of all-electric drive naval vessels that are now being 
planned. 

The Infrared (IR) FEL is a multi-pass, high-power laser operating in the 
near infrared at a nominal wavelength of X= I5(im and micropulse frequency of 
7.5 MHz [1]. The peak current is 72 A with a short micropulse length of 
/ s = 0.05 cm. The electron beam energy is (y-1 )mc 2 = 45 MeV where the 
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Lorentz factor is y=89, the electron mass is m, and c is the speed of light. 
The high-quality electron beam has energy spread of only Ay/y = 0.001 with a 
normalized emittance of e„ = ye = 15 mm-mrad. [11] 

An FEL with a high-quality electron beam operating at IR wavelengths 
tends to have an electron beam that is much smaller than the optical mode in 
the transverse dimension, and smaller than the slippage distance in the 
longitudinal dimension. The optical mode waist radius is proportional to X 1/2 , 
and the slippage distance is proportional to X. In the transverse direction, mode 
distortion tends to increase the gain, while in the longitudinal direction, the 
short-pulse slippage effects tend to decrease the gain. [11] 

In the CEBAF case, the electron beam radius is r e = 0.03 cm, and is much 
smaller than the radius of the optical mode, tv 0 = 0-' | 5 cm. Simple theory 
indicates gain is proportional to the product of filling factor F and dimensionless 
beam current j, but the actual gain observed in simulations is found to be a 
function of filling factor F. [10,11] 

B. CEBAF IR FEL PARAMETERS 

The length of the proposed linearly-polarized undulator is L = 1.5 m with 
wavelength of Xq = 6 cm over N = 25 periods. The linearly-polarized undulator 
peak field strength is B = 4.4 kG, giving an undulator parameter of 
K = eB\ 0 /2'f2Kmc 2 = 1.76. A non-focussed electron beam with smooth 
propagation and minimum electron phase velocity spread is desirable for 
injection into the undulator field. This requires a matching condition between 
the rms radial beam spread, r, and the rms angular spread §, such that 
Kk 0 r = V2y 0 [9]. The emittance is e = re, so that the matched electron beam 
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radius is then r e = (yeX 0 /V2jcK') 1/2 = 0.03 cm. The off-axis injection caused by the 
radial and axial spreads, and the off-axis focusing fields cause a slow, 
transverse "betatron" oscillation of the electrons, with NK/ y=0.5 betatron 
oscillations along the undulator [9]. The amplitude of the electron’s transverse 
wiggling motion is Ax = KXq/^2k y = 0.03 cm, and is comparable to the size of the 
electron beam. The nominal resonant wavelength is X = X 0 (1+K' 2 )/2y 2 = 15 pm, 
but the CEBAF IR FEL will be tunable from 5pm to 20pm. The dimensionless 
electron beam current j = %N[neKUJfprfmc 2 = 138, where p = 4x10 12 cm~ 3 is 
the electron particle density, and the Bessel function factors, 
JJ = Jo(5Wi(5) = 0.78 with % = K 2 /(2(1+K 2 )) = 0.38, describe the reduced coupling 
in a linearly-polarized undulator. [11] 

The value of the Rayleigh length that minimizes the optical mode volume 
is Z 0 = L/VT 2 = 43 cm [9]. Assuming the optical mode is shaped as a lowest- 
order Gaussian with its waist at the center of the undulator {x w = 0.5), the filling 
factor averaged along the undulator’s length is F = r 2 /[tV 0 2 (l+L 2 /12Z,?)] = 0.027 
where the area of the optical mode waist is kW$ = ZqX [9]. The theoretical 
single-mode gain is given by G 0 = 0.135 jF = 50 % [12], and does not include the 
effects of beam quality, self-consistent field growth, short-pulses, or optical 
diffraction. When the interaction includes self-consistent growth of the optical 
field, the gain increases to G = 57 %. The electron beam energy spread has 
the minor effect of reducing the gain by only 0.5 %, and the effect of emittance 
is even smaller. The CEBAF IR FEL electron beam is of such high quality, we 
will not include the effects in further simulations. [11] 
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C. SHORT-PULSE EFFECTS 



When the length of the electron pulse, l g , is comparable to the slippage 
distance, As = NX, the Fourier components of the pulse current are comparable 
to the gain bandwidth, and short-pulse effects influence the interactions in the 
undulator [9]. As an optical pulse bounces off the resonator mirrors and arrives 
at the entrance to the undulator, at t = 0, an electron micropulse from the 
accelerator is timed to arrive simultaneously. A slight displacement of one of 
the resonator mirrors, termed the "desynchronism” d = As/nX, alters the 
distance the rebounding light must travel, and therefore will serve as a means 
to adjust the arrival of the optical pulse relative to the subsequent electron 
pulses. Desynchronism, d, is normalized to the slippage distance, and an exact 
synchronism between arrivals of the light and the electron pulses occurs when 
d = 0. 

For the CEBAF IR FEL operating at the nominal wavelength X=l5^m, 
Figure 4-2 shows an undulator with exact synchronism d= 0. The window for 
calculation, at middle left, is four slippage distances long, and travels with the 
light pulse at speed c. The electron pulse, traveling at slower speed (5 z c, slips 
backward through one-fourth of the window per pass. The positions of the 
short electron pulse relative to the optical pulse is shown at left bottom. The 
electron pulse position y'(z^c), initially is shown as the lighter parabolic shape at 
x = 0. As bunching develops the electron pulse slips backward, and the final 
position is indicated by the darker shape at x = 1. The pulse form is parabolic 
shape y'(z) = y(l-2z 2 /of) for y'(z)> 0, or otherwise zero. The electron normalized 
pulse length o 2 = l e /NX = 1.4, and the peak current is \F = 3.8 where F = 0.027. 
The optical amplitude grows proportionally to j on each pass, while amplitude 
losses due to the resonator are calculated by e~ vo for each pass. [9,1 1] 
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(d = 0 ). 
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In order to start the optical pulse growth, random shot noise is generated 
by addition of a small phase, 5£ = 0.0001 , to each of the sample electron 
phases. The optical pulse is free to evolve subject to gain effects of the 
electron pulse, and resonator losses due to Q over n = 800 passes through the 
undulator. Each pass, a new electron enters the undulator to drive the 
rebounding optical pulse. The shot noise creates a low power, broadband 
optical pulse in the first few passes with a length that is roughly equal to the 
sum of the slippage distance and the electron pulse length, c 2 + 1. This makes 
sense because this is the end-to-end distance over which the optical pulse is 
exposed to the amplifying electrons when at perfect synchronism. 

The optical field is amplified at the rear of the optical pulse centroid after 
n = 10 passes. The electrons are amplifying the optical pulse, which has 
partially passed over the electrons. It is clear that as the number of passes, n, 
increases, the centroid of the optical pulse moves backward (to the left) due to 
the continual amplification at the rear of the pulse. The amplitude evolution 
peaks near n = 200 passes, and then decreases. The final shape of the pulse, 
la(z,n)l, is shown, top-left, above the amplitude evolution plot. The maximum 
amplitude of the optical field is \a{z,n)\ = 58 is shown on the grey scale in 
black. A single contour line signifies half-maximum values, while zero 
amplitude is in white. 

The power spectrum evolution (middle-center) starts near peak gain at 
v 0 = 2.6, and grows in width at n = 150 passes. The broadened range of phase 
velocities v at this point is caused by a short pulse. The weak-field gain 
spectrum, G(v), is plotted on the same scale below for reference. The final 
power spectrum, P(v,n), shows the spread in phase velocities for the last pass, 
n = 800. The square tick-mark at the top of the power spectrum indicates the 
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center of the final power spectrum, while the pointed tick-mark shows the 
central wavelength of the initial radiation at resonance. The power evolution, 
P(n), also shows rapid increase, peaking near n = 150, then tapers ultimately to 
zero power. The light wave in exact synchronism systems eventually drifts 
away from the slower electrons leaving the FEL with steady-state power of 
zero. 

The electron beam phase velocity evolution (right-center) shows a single 
phase velocity emerging because of the vanishing optical power. On the final 
phase velocity plot, f(v,n) (right-top), the coincident initial and final phase 
velocity tick-marks indicate near-constant energy of the electron pulse due to a 
optical field. 

In an undulator with desynchronism d > 0, the centroid of the optical pulse 
arrives at i=0 slightly ahead of the centroid of the electron pulse which then 
passes back over the optical pulse. As the pulses propagate and x increases, 
the electrons bunch and radiate, amplifying the rear portions of the optical pulse 
with gain. The shape of the optical pulse is distorted with each successive 
pass through the undulator. Because the light pulse essentially grows from the 
rear, the light pulse will be distorted in amplitude over the slippage distance as 
the rear of the pulse grows. Figure 4-3 shows this effect with a small 
desynchronism d = 0.005, for the CEBAF IR FEL described in Fig 4-2. 

In Fig 4-3, altering the FEL with small desynchronism d = 0.005 results in 
major changes is performance of the laser system. The optical amplitude 
evolution survives and produces an amplitude in steady state at n = 300 passes 
which has the dual-peak modulation in the final pass window la(z,n)l. 
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**** pel Pulse Evolution **** 
jF=3 . 8 a z = 1.4 d=0 . 005 8^=0.0001 

0=20 N=25 




-2 z 2-19 v 19 0 n 800 



Figure 4-3: An FEL system with desynchronism and amplitude 
modulation due to short-pulse effects. 
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The power spectrum evolution, P(v,n), and optical power, P{n), also 
achieve steady state, with the maximum power much greater for this laser with 
non-zero desynchronism than the earlier synchronous FEL. 

Though the final optical pulse shape \a{z,n)\ is frozen in a single picture, 
the series of pictures for each successive pass shows the broad amplitude 
subpulses originate in an area of z~- 1, and grow while moving to the right. 
This effectively describes the optical pulse being amplified along its rear, and a 
subpulse passing forward away from the slower electrons. The magnitudes of 
the subpulses would depend on the parameters of the system under simulation. 
The amplification and growth of the rear of the optical pulse also has the result 
that the centroid of the pulse moves at a speed less than c, which effectively 
lowers the pulse group velocity. 

The tick marks in the final power spectrum picture shows the spread of 
center phase velocities expected, and there is a small sideband caused by 
trapped-particle instability. The wider spread of the electron spectrum f(v,n) is 
proportional to the separatrix peak-to-peak height, and the initial and final tick 
marks show a drop in electron phase velocities, indicating a net transfer of 
energy to the optical field. 

The optical wavelength of an FEL is defined X = Ao((l + K z )i2-f). In order 
to tune an FEL to a new wavelength, the undulator parameter, K = eEX 0 /2nmc 2 , 
is altered by changing the mean magnetic field E. This in turn changes the 
optical mode waist, W 0 , the filling factor, F, the dimensionless electron beam 
current, j, the dimensionless pulse length, o z = l e /NX, and several other FEL 
parameters. The effects short wavelengths for the CEBAF IR FEL operating at 
X = 5|am is shown in Fig 4-4. 
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*+** FEL Pulse Evolution **+* 



j=11.5a z =4.3 d=0 . 005 5^=0.0001 
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Figure 4-4: An FEL operating at X = 5^m with long electron pulse 
length a z . 
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The dimensionless pulse length here is a z = 4.3, three times the length of 
the electron pulse for the FEL operating at X=15^m shown in Fig 4-3. The 
amplitude evolution plot clearly shows amplitude spikes growing from the rear of 
the optical pulse, and moving forward over several passes. The longer pulse 
results in a longer optical pulse spread along eight or ten slippage distances 
and la(z,n)l amplitude spikes spaced roughly one slippage distance apart. In 
contrast to the narrow power spectrum of the short pulse FEL of Fig 4-3, here 
the power spectrum is spread among several sidebands without a steady-state. 

Power evolution, P(n), rises quickly but never achieves a steady value. 
The higher average power is a direct result of the longer electron pulse 
interacting with the optical pulse for a longer period and so transferring more 
energy. The electron spectrum, f(v,n), of Fig 4-4 is again broad, with the 
average phase velocity dropping due to net energy transfer from the electrons 
to the optical pulse. 

While electron pulse described by a z = 4.3, shown in Fig 4-4 has the 
smallest slippage within the tuning range 5|xm s X < 20|am, the largest slippage 
is o 2 = l.l at 20nm. This system is shown in Fig 4-5. Figures 4-4 and 4-5 
serve well to contrast the various effects of pulse length on system 
performance. Short-pulse effects generally result in short, narrow optical pulses 
spread over fewer slippage distances and comprised of fewer spikes. 

Figure 4-6 shows an FEL operating near X=15|im with relatively high 
resonator losses, 0 = 50, and a moderate desynchronism d = 0.0165. This 
higher value of d causes the modulation of the amplitude evolution to stand out 
clearly. [1 1] 
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**** FEL Pulse Evolution **** 
jF=2 . 9 a 2 =l.l d=0 .005 6^=0.0001 
0=20 N=25 




-3 z 3 -19 v 19 0 n 800 



Figure 4-5: FEL operating at X = 20n.m and short electron pulse length 
o z = 1.1. 
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**** pex. Pulse Evolution **** 
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Figure 4-6: An FEL with steady state limit cycle behavior in the power 
evolution P{n). [1 1] 
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The power evolution P(n) exhibits sharp overshoot followed by a steady 
state oscillation in power as new optical pulses grow periodically. The period of 
the oscillations appears to be An = 90 passes. This oscillation is known as 
"limit-cycle behavior." This phenomena occurs in regimes of moderate 
desynchronism where a stable overlap does not occur between the optical 
pulse and the electron pulse. Moderate field strengths cause about one 
synchrotron oscillation, which causes the optical pulse to continually shift from 
one shape to another. [9,1 1] 



D. CEBAF IR FEL TUNING RANGE AND GAIN 

The tuning range of the CEBAF IR FEL (5|j.m £ X£20|j.m), combined with 
the effects of varying desynchronism on gain, allows for a broad range of 
system performance. The broad tuning range of the FEL includes electron 
pulse lengths ranging 1.1<o 7 <4.3, which would imply a wide range of system 
gain. Neglecting the effects of the slight energy spread, and using multiple 
passes to achieve steady-state gain, simulations were run over a range of 
desynchronisms. 

The shortest wavelength of the tuning range, X = 5^m, has the longest 
pulse length c 2 = 4.3. Varying desynchronism d from 0.005 through 0.06 
resulted in a steady single-pass gain of G = 170 %. For wavelength X = 10n.m, 
with pulse length o z = 2.2, the gain is only about 60 %. At X = 15|im, the gain is 
about 35 %. At the longest wavelength of this tuning range, X = 20jam, the gain 
is about 20 %. Fig 4-7 shows the variation of the system gain across the 
wavelength tuning band of the CEBAF IR FEL. 
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Figure 4-7: Gain as a function of wavelength \ for the CEBAF IRFEL. 



The longest wavelength A. = 20^m has the lowest gain in the tuning 
bandwidth with G = 22 %. We considered lengthening the undulator from 
N = 25 to N = 35 periods in order to increase gain. The dimensionless current j 
increased from 138 to 379. The undulator length L increased from 150 cm to 
210 cm, while the dimensionless electron pulse length decreased from a z = 1.1 
to 0.77. The average filling factor F was reduced by half, but since j increased, 
the product jF increased by a factor of three. While increasing jF increases 
gain, the decrease in o z decreases gain. The net result from the lengthening of 
the undulator is determined by simulation. The gain is improved to G = 36 % 
for N = 35 periods and G = 50 % for N = 45 periods. The cost of increasing the 
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length of the undulator in order to improve gain at longer wavelengths would not 
be large. 
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V. OPTICAL MODE DISTORTION 



A. BACKGROUND 

After originating from spontaneous emission, the optical field in the 
undulator is a sum of rebounding light waves superimposed with newly 
generated light from the bunched electrons within the electron beam [15]. The 
longitudinal effects of short pulses were described in Chapter IV, and the 
transverse effects caused by natural diffraction and distortion by the electron 
beam will be treated in this chapter. The effects of short pulses will not be 
included in this chapter. For the CEBAF IR FEL F = 0.027 [10]. As electron 
bunching develops late in the length of the undulator, the new light is created in 
an area that is = F times the size of the mode area. The newly created light 
from the small electron beam, has a short Rayleigh length and diffracts rapidly. 



B. GENERAL MODE DISTORTION 



The effects of the small electron beam size acting on the optical mode 
may be simulated through numerical integration of the parabolic wave equation, 
3-9, coupled with the Lorentz force equation, 3-1. To calculate the effects of 
mode distortion, it is necessary to shift the transverse dimensions into a 
dimensionless coordinate system, normalized to the mode size, x = xV/c/2 L . 
The optical amplitude, la I, and the dimensionless current density, j , are also 
shifted into this coordinate system. The parabolic wave equation now becomes 




\ 

dz 



a(xVc) = -</e' ; >, 



5-1 
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where the operator V 2 = d 2 /dx 2 + d 2 /dy 2 , properly describes the diffraction of the 
optical wave in transverse dimensions [13]. An important effect of diffraction is 
self-focusing of the optical mode back into the electron beam. The optical 
phase shift, <|>(x), occurs in the center of the optical mode, near the electron 
beam, and acts to focus the light back onto the undulator axis. The FEL 
interaction causes the light to remain closer to the electron beam than would 
occur in the presence of natural diffraction alone. This effect is greater when 
gain is larger. When the optical mode is larger than V2L/A the distortion 
operator = 1, and can be neglected. [13] 



C. MODE DISTORTION SIMULATIONS 

Figure 5-1 shows a numerical simulation of a Gaussian optical wavefront 
propagating without the effects of electron beam current (j = 0), as x goes from 
0 to 1 in an undulator with N = 25 periods. The wavefront has weak initial 
optical field amplitude, a 0 = l.O, at mode center, and mode curvature giving 
normalized Rayleigh length of z 0 = 0.5 focussed at \ w = 0.5. Transverse 
dimensions are normalized to the mode waist <k/2L . The grey scale at upper 
right shows the intensity for minimum (white), and for maximum (black), with 
five white contour lines to illustrate the focusing and diffraction of the light wave. 

The optical amplitude evolution plot, la(x,x)l, (upper-left) shows the optical 
mode evolution over a single pass through the undulator. The five contour lines 
show the focusing and diffraction expected for an undistorted optical wavefront. 
The optical amplitude is greatest at the mode center, and decreases radially as 
shown in the transverse optical amplitude plot, la(x,y)l, representing the final 
mode at the end of undulator. The optical phase shift, <t>(0,x), shifts negatively, 
A<(> = -x/z 0 , during diffraction. 
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Figure 5-2 shows an FEL with the same optical parameters as shown in 
Fig 5-1, but with optical mode distortion caused by an electron beam current 
j = 138 for the CEBAF IR FEL The electron beam has normalized radius, 
a x = r b /'lLX/K = 0.125, and has initial phase velocity v 0 = 4.6. With a single 
contour line, the optical amplitude plot, la(x,x)l, shows focusing of the light with 
Rayleigh length z 0 = 0.5 in toward the mode waist positioned at x w = 0.5. 
Relative to the undistorted field of Fig 5-1 , the field distortion in Fig 5-2 is clear. 
The maximum amplitude, la(x,x)l, is greater, and the contour line spreads 
rapidly from the beginning of the undulator as the distortion of the electron 
beam drives up the amplitude of the optical mode along the second half of the 
undulator. The final transverse amplitude plot, la(x,yl, shows the relatively 
narrow mode focussed around the electron beam. The electron bunching 
current, o(x,x) « j(r)< cos(C, <t>)>, evolves as the electron phases, C. bunch within 
each optical wavelength. The plot of final electron bunching, o(x,y), shows the 
relative disparity in size of the electron beam relative to the final optical mode, 
la(x,y)l. [11] 

The electron phase velocity evolution, f{v,x) in Fig 5-2, starts at v 0 = 4.6, 
and broadens slightly due to bunching. There is also slight bunching shown in 
the phase-space plot, (C,v). The optical phase shift, <t>(0,x), is also affected by 
distortion effects of the beam current, From 5-1, the FEL optical phase shift 
1 = y<sin(C + 4>)>/la I, is maximized at (C + <t>) = nJ2. The electron phase, C. is a 
function of dimensionless time, x, so that for strong-field FEL interaction, 
<t> = y'x/la I, and for weak fields, <t> = yx 3 /1 2. [9] 
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Figure 5-1 : Non-distorted optical mode propagation diffraction. 

In either a strong or weak fields, the FEL interaction gives a positive shift 
to the optical phase. Variation of <t>(0,t) as i increases from 0 to 1 is driven by 
electron beam interaction. 

Exponential growth of gain and power is delayed in the undulator until the 
bunching time, i g = {2/j) w = 0.25, has passed, after which the strong beam 
current, j, drives up the gain and power [9], The gain for the FEL in Fig 5-2 is 
ln(l + G(-c)) = 0.846 or G = 133%. Theoretical single-mode gain for this system 
predicts only G 0 = 0.135/F = 43%. In single-mode theory, the dimensionless 
current density, /, is proportional to the electron beam density, p, and is 
therefore inversely proportional to the square of the electron beam radius, or 



;«r 9 - 2 . [11] 
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Figure 5-2: FEL with optical mode distortion. [11] 
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The filling factor, F, is proportional to the electron beam radius, F « r e 2 , 
and so the theoretical gain, G 0 « jF, becomes independent of the electron 
beam size as long as the electron beam is contained within the optical mode. 
But single-mode theory does not account for the mode distortion caused by the 
intense current and small filling factor of the FEL in Fig 5-2. [11] 

Figure 5-3 shows the effects of shifting mode waist position, i w , and 
Rayleigh length, z 0 , to achieve increased gain. Showing strong diffraction and 
distortion effects, the amplitude evolution, la(x,x)l, shows the light focussed 
inward to x w = 0.19, and then sharply diffracting outward from a narrow waist 
near x = 0.40. This is due in part to the short Rayleigh length, z 0 = 0.14, which 
causes light to diffract rapidly. The plot of the final mode, la(x,y)l, clearly 
shows the narrow focusing of the optical mode around the electron beam, and 
the rapid spread of highly diffracted light. 

The larger, irregular shift of <j>(0,x) for the FEL in Fig 5-3 is due to the short 
Rayleigh length, z 0 = 0.14, and the varying effects of the field amplitude la I 
along the undulator. The apparent discontinuity late in the plot of 4>(0,x) 
represents the optical phase shift decreasing to = -n. Electron bunching 
occurs later in the undulator than for the FEL of Fig 5-2, as shown by the 
bunching current plot, o(x,x). The final gain is much greater than for Fig 5-2, 
with ln(1 + G(x)) = 1.14, or G(x) = 2l2%. Single-mode theory, neglecting the 
effects of mode distortion, predicts this FEL to have gain G 0 = 40 %. 
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Figure 5-3: An FEL with short Rayleigh length and high distortion. 
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D. MODE PARAMETERS, RESONANCE AND GAIN 

Variation of Rayleigh length and mode waist position and the resultant 
gain curves are shown in Figures 5-4 and 5-5. Typical values for Rayleigh 
length and mode waist position are z 0 = % Wq/LX = Zq/L = I/V 12 or z 0 = 0.5, and 
x w = 0.5, but simulations for the CEBAF IR FEL produced peak gain at 
z 0 = 0.33. Using this Rayleigh length, the peak gain for variation of mode waist 
position was located at x w = 0.35. 

Figure 5-6 shows the optical gain spectrum, G(v 0 ), plotting the gain 
against the FEL phase velocity, v 0 = L[(k+k 0 )$ z -k] [7], Rayleigh length, 
z 0 = 0.14, and mode waist position, t iv = 0.19, are selected to maximize the 
single-pass gain for the small CEBAF electron beam. When X ~ 15 urn, the 
CEBAF IR FEL is near resonance v 0 = 0, and Fig 5-6 shows the gain around 
resonance. When the phase velocity is off-resonance by Av 0 = re, the shift in 
wavelength is AX/X ~ A\qI2kN -V2N-2 %. The electron beam has a radial 
parabolic shape described by j{r) = ;'(1-r 2 /2af) for r < V 2 a g , and zero otherwise. 
The peak current is j = 138, and the normalized beam radius is 
o e = rg'fjUIJ. = 0.125. [11] 

The simulation in Fig 5-6 calculates gain including the transverse effects 
of diffraction, but ignores the longitudinal short-pulse effects. At each value of 
v 0 , a wavefront enters the FEL undulator with peak field strength a 0 = 1 in the 
weak field regime, where a 0 = 4nNeKLEi-fmc 2 , and E is the optical electric field 
strength. The peak gain in Fig 5-6 is G = 233 % at phase velocity v 0 = 5. The 
region of positive gain is larger than the region of negative gain, as is typical of 
high current FELs [9]. The gain spectrum tails off smoothly as v 0 increases on 
the right of Fig 5-6. [11] 
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Figure 5-4: Variation of gain with Rayleigh length, z 0 . 




Figure 5-5: Variation of gain with mode waist position, z w . 
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*** FEL 3d simulation, single-pass gain *** 
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Figure 5-6: The optimized CEBAF IF FEL gain spectrum including only 
transverse modes. 
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E. GAIN AND SMALL ELECTRON BEAMS 



The FEL gain observed for Figs 5-2 and 5-3 indicates that single-mode 
theory can significantly underestimate the effects of diffraction and distortion on 
actual FEL gain. Calculating the single-pass gain over a broad range of filling 
factors, the dependence of FEL gain on electron beam size is examined in Fig 
5-7. Neglecting optical mode distortion, the single-mode FEL gain is 
G 0 = 0.l35y'F. The gain ratio G/G 0 is determined numerically and plotted as a 
function of the electron beam size compared to the optical mode waist, a e /w 0 , 
with a single line on the plot signifying G/G 0 =1.0. The product of the peak 
beam current and the beam area, yof ~ jF, is held constant at yF = 0.l. The 
electron beam shape is taken to be parabolic, j{r) = j{\-r 2 t2ol) for r < <2 o e and 
zero otherwise. A small value of jF is used to insure that there are no high 
gain effects for values on the left side of Fig 5-7 where F = oJ/wq is small and y 
may be large. The Rayleigh length is, z 0 = 1/VT2 = 0.29, for optical single-mode 
gain G 0 and position of the mode waist is, x w = 0.5. These values are chosen 
to describe a more typical FEL mode instead of the optimized CEBAF IR FEL. 
The phase velocity is v 0 = 5 for peak single-mode gain [12], and the peak 
optical amplitude is a 0 = 1 for weak fields. [11] 

The plot shows a clear dependence on the electron beam size at all 
values plotted. The single-mode theory is correct for the range of values o e /w 0 
where the G = G 0 . We see that there is no range of values where the single- 
mode theory is correct, but only a single point. For a small electron beam, the 
single-pass gain including many modes can be significantly higher than the 
single-mode gain, G 0 . It is emphasized that the mode distortion crucial to the 
increased value of G/G 0 is not caused by high current. [11] 
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*** FEL 3d simulation, single-pass gain *** 
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Figure 5-7: General gain curve for FELs with small electron beams. 
[ 11 ] 
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When the product jF is decreased by xlOO to jF = 0.001, the plot in Fig 5- 
7 remains substantially the same. When jF is increased to unity, or above, the 
gain enhancement due to mode distortion is further increased by high-gain 
effects on the left side of Fig 5-7 where beam radius is small and j is large. It 
is not clear from this analysis what part of the distorted mode will remain in a 
resonator after many passes. The actual gain of an FEL oscillator may be 
reduced from that shown in Fig 5-7 depending on the resonator design 
[13,15,16]. No known work explains Fig 5-7, and in particular, its independence 
of the value of jF below unity, where there is vanishingly small mode distortion. 
[ 11 ] 

F. PROPAGATION AND INTERFERENCE 

A modification of the simulation that produced Figures 5-1 to 5-3 subtracts 
the input optical field from the total field to allow analysis of the propagation of 
newly created light during a single pass through the undulator. This process is 
interesting because of the complex optical field interactions shown in Fig 5-8. 
Because the amplitude plots, Ala(x,x)l and Ala(x,y)l, are generated by 
subtraction, it is possible to have negative amplitude values, representing 
destructive interference in the total optical field, reducing field amplitude. 
Destructive interference occurs when "new" light from the bunching of electrons 
encounters "old" light with phase difference of Ia<(>I = Ax/z 0 = tc radians. Noting 
that the newly generated light appears in the amplitude evolution, Ala(x,x)l, near 
Ax = 1/3, and holding optical mode waist, w 0 = [z 0 + 1/12z 0 ] 1/2 , constant by the 
input Rayleigh length, z 0 = w£ = o|, a phase shift relationship may be found: 
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Figure 5-8: Negative interference of propagated light. 
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5-2 




For electron beam radius, o e £ 0.3, conditions for destructive interference may 
exist. 

The plot of propagated optical amplitude evolution, Ala(x,x)l, shows an 
area of destructive interference starting near x = 0.7, which is obvious in the 
final amplitude plot, Ala(x,y)l, as a light area in the center of the optical mode. 
To evaluate the interference effect over a range of values of y'oj, the product 
was held constant while y and o e were varied. Destructive interference was 
observed for 0.12 <c s < 0.3, and the point of initial destructive interference 
occurred earlier in the undulator as the electron beam, o e , is narrowed. This is 
in part caused by the constant product y'o|, which rapidly raises beam current, 
y, as the electron beam narrows. This appears as a radial ’dimple’ across the 
propagated wavefront after a single pass through the undulator. In order to 
determine whether the phenomenon of destructive interference in the optical 
wavefront continues outside the undulator, a multi-pass simulation within the 
entire resonator is required. This destructive interference over a single pass 
was first investigated by LT Carl Bice in his thesis in December, 1991, and is a 
continuation of that work.[14] 
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VI. MULTI-PASS OPTICAL MODE DISTORTION 



A. BACKGROUND 

To continue the investigation of mode distortion, the single-pass three 
dimensional simulation used in Chapter V is modified to calculate the transverse 
beam effects over multiple passes. This enables further research into the 
interference effects shown in Fig 5-8, and the relationship between gain and 
electron beam size as shown in Fig 5-7. The multiple pass simulations 
introduce optical mode selection caused by the finite size of the resonator 
mirrors. Through many passes, highly diffracted modes may be lost due to 
mirror size and placement. 

B. MULTIPLE PASS PROPAGATION AND INTERFERENCE 

Investigation of the single-pass optical interference pattern shown in Fig 
5-8 over multiple passes will be affected by the number of passes, n, and the 
mirror design. Figure 6-1 shows the three dimensional simulation modified for 
multiple passes. For this regime, the initial optical amplitude, a 0 = 0.05, for a 
weak field system. At either end of the amplitude plot, la(x,n)l, are curved 
mirrors, spaced x m = 2 apart, or twice the undulator length. The separation 
distance of a typical FEL may be larger, but this choice allows ease of 
computation. The dimensionless mirror radius of curvature, r c , is determined by 
the mirror spacing, x m , and the normalized Rayleigh length: r c = x m /2 + 2z$/x m . 
The actual mirror radius of curvature is given by Lr c . The mirror losses are 
described by Q = 50, and the mirror edge-loss is k = 0.001. 
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FEL WAVEFRONTS (HOLE-COUPLED OUTPUT) 
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Figure 6-1 : Three dimensional simulation for multiple passes and 
destructive interference. 
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Edge-losses occur when the mirror radius, r m , is not wide enough to capture 
the higher-order optical modes, losing them to the non-reflecting sides of the 
resonator. The dimensionless mirror radius is a function of the Rayleigh length, 
mirror spacing and edge loss, r m = [z 0 (1 + x£/4z^) (1/2 log((1 + k)/k))] 1/2 . The 
amplitude plot, la(x,n)l (upper-left), shows the evolution and growth of the 
optical mode through n = 16 passes, as scaled along the base of the plot. The 
single contour line shows that the amplitude begins to grow in the last few 
passes. The plot of resonator mode evolution, la(x,x)l (center), shows the 
entire optical mode at pass n = 16. The undulator ends are represented by tick 
marks at x = 0 and x = 1 along the plot base. 

Careful analysis reveals that a radial "dimple”, similar to that found in the 
single-pass interference plot of Fig 5-8, is present near the right end of the 
optical mode. In the grey scale used, the destructive interference described in 
Chapter V is present as two unshaded rays that propagate from just outside the 
undulator markings toward the right-side mirror. The final transverse amplitude 
plot, la(x,y)l, shows the amplitude at the last pass, n = 16 at the end of the 
undulator. 

The electron phase velocity evolution, f(v,n), shows only small modulation 
in weak fields. The phase-space plot, (£.v), shows slight bunching in this weak 
field regime. Power evolution, P(n), starts to grow at approximately the same 
point as the amplitude in la(x,n)l. Gain, G, reaches a steady-state near n = 8, 
where G = 55%, and single-mode gain.G 0 « jF predicts only 28%. 

Figure 6-2 shows an FEL similar to that of Fig 6-1 , but in the strong field 
regime, starting with a 0 = 5.0, and running over n = 160 passes. The amplitude 
evolution plot, la(x,n)l, has a much greater maximum of 39. 



64 



FEL WAVEFRONTS (HOLE-COUPLED OUTPUT) 
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Figure 6-2: Three dimensional simulation without destructive 

interference. 
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The final transverse amplitude plot, la(x,y)l, and the resonator plot, 
la(x,x)l, show an optical mode with no radial ’’dimple". The phase evolution 
plot, f(y,n), and phase-space plot, (£,v), show complex bunching due to strong 
field saturation. Power evolution, P(n), starts near where the amplitude plot 
begins to grow, and continues to increase through the n = 160 passes. If this 
simulation continued for more passes, the power evolution, P(n), would achieve 
a steady state value. The gain plot shows an early peak, followed by a quick 
reduction to near-zero steady-state near n = 50. It appears that the interference 
"dimple" is not present at saturation in the FEL oscillator. 

C. GAIN AND SMALL ELECTRON BEAMS OVER MULTIPLE PASSES 

To investigate the effects of electron beam size on gain over multiple 
passes, the product y'of is held constant, while the electron beam radius is 
decreased. Actual FEL gain, G, was numerically calculated over several 
passes, while the single-mode gain is given by G 0 = 0.135 jF « /oj. In order to 
determine whether the single-pass relationship of gain on electron beam size 
shown in Fig 5-7 applies for multiple passes through the resonator, simulations 
were run to evaluate gain. For these simulations, mirror loss are Q = 1000, and 
the edge loss is k = 0.001. The Rayleigh length, z 0 , edge loss, k, and mirror 
spacing, x m , result in mirror radius, r m = 2.7, and radius of curvature r c = 0.9. 

Figure 6-3 shows the ratio of the actual gain to the single-mode gain, 
G/G 0 , plotted against the ratio of dimensionless beam radius to mode waist, 
o e /w 0 while holding yo| = 0.4. Fig 6-4 shows a similar plot for the value 
jog = 0.6. To prevent the influence of high-gain effects, the beam radius is 
limited by o e /w 0 > 0.1. 
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Figure 6-3: Multiple-pass gain curve for /of = 0.4. 




Figure 6-4: Multiple-pass gain curve for /of = 0.6. 
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The curves resemble the single-pass results in Fig 5-7, with the gain 
ratios decreasing steadily as the electron beam radii widen toward the size of 
the mode waist. The single-mode gain product yof is greater for Fig 6-4, so the 
gain curve for Fig 6-4 is higher than for Fig 6-3. Similar to Fig 5-7, there is no 
range of values where single-mode theory is correct, only single points where 
G/G 0 = 1 . 
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VII. CONCLUSIONS 



The development of the FEL as a useful tool continues through the 
progress of theoretical and experimental efforts in several countries. The 
potential for military application of the FEL are more promising than other 
sources of directed energy. A major challenge to all high energy light sources 
is the complex effects of atmospheric propagation. Innovative beam relay 
methods may enable these obstacles to be minimized. 

The proposed design of the CEBAF IR FEL will make important 
contributions toward realizing the goal of a high power, high efficiency source of 
coherent radiation for applications in medicine, industry, and possibly the 
military. Simulations of the CEBAF IR FEL in Chapters IV and V have 
contributed to the FEL design selection process. The broad tuning range of the 
CEBAF IR FEL was analyzed for system gain and power, and the influence of 
short electron pulse effects. The small electron beam of this system created 
strong distortion effects which were simulated in Chapter V. The interaction of 
the optical mode and the intense beam current resulted in higher FEL gain than 
predicted by single-mode gain theory. There exists a complex destructive 
interference in the optical mode that results in a radial "dimple" in the 
propagated wavefront. This phenomenon, and system gain, was explored in 
detail for single pass simulations in Chapter V and multiple passes in Chapter 
VI. 
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